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Abstract 



> 

^-j. | The structure of the collective low-energy J T = CP (T=0 and T=l) modes is 

o ' 

CO | studied for a doubly magic nucleus in a schematic analytic model of RPA. The 



phonon states (T = 0, 1) lie at energies Et=o(0 ) < u and Et=i(0 ) > CJ, 
where oj is the oscillator frequency. The matrix elements of P-odd and P- 



and T-odd weak one-body potentials connecting the ground state to these 

CT-states, W co u, are enhanced by the factor ~ 2(^) 1 / 2 A 1 / 3 ~ 10 as compared 

^ ; to the single-particle value w sp what can result in values | W co u | ~ 20 — 30eV 

if standard values of DDH parameters are used for w sp . Similar enhancement 

arises in the P- and T-odd case. 
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The purpose of this work is to study low-lying collective nuclear states with quantum 
numbers J n = 0~ and the matrix elements of parity nonconserving potentials connecting 
these states to the ground state. The work is motivated not only by the interest in the 
low-energy 0~ modes, which already attracted attention earlier in relation to pion (pre- 
) condensation, but also by the possibility of the collective enhancement of P- and P-,T-odd 
matrix elements of the pseudoscalar "weak potentials" in nuclei. The importance of collective 
CT states for these effects has been already revealed within doorway state approach to the 
P- violation in compound nuclei Jl| , [g] . It was pointed out that the 0~ component of the 
spin-dipole state is an important mediator of parity mixing between nuclear states. This 
applies to both parity violating (P-odd) and time reversal conserving (T-even) as well as to 
P-odd and time reversal violating (T-odd) interactions. The one-body potentials for these 
two symmetry violating interactions can be approximately described by the operators {Bp) 
for P-odd T-even and (ar) for P-odd, T-odd. When acting on a nuclear state both operators 
will excite the 0~ mode. Since the weak interaction does not conserve isospin one considers 
both isoscalar (T = 0) and (T = 1) J = 0~ spin-dipole states. 

In Refs. renormalization of the weak P- and P-, T-odd potentials due to the residual 
strong interaction were calculated and the positions of low energy 0~ modes was discussed. 
The properties of 0~ (T = 1) states were studied earlier in an elaborated Hartree-Fock+RPA 
approach in Ref. 0. Spin-dipole modes were studied in 0, || within an approximate RPA 
approach. An intersting analysis of related effects is presented in Ref. 0. 

Here we show, using an approximate analytical RPA calculation, that a realistic residual 
nucleon-nucleon interaction which has been widely used to describe various nuclear prop- 
erties does indeed lead to formation of such collective 0~ modes. We give approximate 
quantitative predictions for their characteristics. We confine ourselve here to the case of 
doubly magic nucleus. The density of states is relatively low. We neglect in our study the 
coupling of the one-particle one-hole 0~ states to all other configurations. 

We write the nuclear Hamiltonian in the form 



2 



H = H + V, (1) 

where V is the residual strong interaction and H is the singe-particle shell-model Hamilto- 
nian of nucleons with momenta p a and mass m moving in a spherically symmetric average 
potential U(r) : 

Ho = J2(^ + U(r a )) 7 U(r) = -\U \ + ^ (2) 



a 



(we use the harmonic oscillator well of the depth |£7o| and neglect for simplisity the shell 
splitting). 

We employ the Landau-Migdal parametrization |10| of the residual strong interaction 
chosing V to be the sum of velocity independent V and velocity dependent V\ components, 
V = Vq + V\. The velocity independent term in the residual interaction is given by 

V^l, 2) = C8{n - f 2 )[/ Q + f Tif 2 + #<?i<7 2 + g'r\T%BxB^ (3) 

where C = 300 MeV ■ fm 3 is the universal Migdal constant. (We use the convention 

The spin- and momentum-dependent component of the interaction V\ will be used in two 
forms: the first one can be written fl2| in the form of a Landau-Migdal spin- and velocity 



dependent interaction, 

Vi(l,2) = ^Cp F 2 {g 1 + fi'ir 1 r 2 )(o ! 1 a 2 ) [pip 2 <Kn _ ^2) + PiS(n - f 2 )p 2 + 

p 2 S (fx - r 2 )px + 5(rx- r 2 )pxp 2 ] , (4) 

Pp is the Fermi momentum. In Eqs.([|fD, the parameters / and g are the strength con- 
stants of order of unity (their values will be given below). Generally speaking, they contain 
r-dependence which is usually parametrized via the density dependence. However, this de- 
pendence is considerable only for the constants / and /' which describe spin-independent 
components of the interaction. In what follows, only spin-dependent components of the in- 
teraction are essential, and the corresponding strengths g[ are density independent. This 
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phenomenological interaction can be viewed as a model of more realistic but also more com- 
plicated forces. It is therefore convenient to use it first because of simplicity of calculations. 
In fact, the results obtained for this interaction manifest the correct tendencies (see also Q). 

We also make calculations for the 7T- and p- meson exchange interaction in its explicit 
form 



k +p = ^(m) 

mi 



f, 



2 



(V x ai)(V x a 2 ) 



r 



(5) 



using it instead of V\ in Eq.(£|). Here, m n (m p ) are the 7r- meson (p-meson) mass, f n and f p 
are the coupling constants, and V is the derivative with respect to the nucleon separation r = 
r*i — r*2; x denotes the external vector product. The standard "Lorentz-Lorentz" correction 
(see, e.g., 



[L3|]) is assumed [|]. 

In the general case, the RPA equations can be written in the form of equations of motion 
for the phonon creation operators A 1 

- E n Ai = [Al n RPA ] = [Al H ] + ([Al V}), (6) 

where the first equality comes from the commutator of the phonon creation operator A^ n 
with the RPA Hamiltonian 1~Lrpa = const + J2 n EnA^An which is derived from the initial 
Hamiltonian Hq + V when the RPA phonons are used. The expectation value in (|]) is 
evaluated using the uncorrelated ground state |0 + ). Within the RPA, the phonon operators 



are 



— ^mi^-wpi ~ ^2 ^mf-i^-m, (7) 
im im 

where the c\ and q, are the operators of the creation and annihilation of nucleons in the 
single-particle states marked by a and b, respectively. The amplitudes X and Y are related 
to the "forward" and "backward" going graphs, respectively. The index n marks the modes, 
the total number of which A/", coincides with the total number of the particle-hole pairs with 
quantum numbers 0~. 
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The 0~ particle-hole states have energies la;, 3a;, ... . We are interested in low-energy 
collective states which should be dominated by the lu transitions. The corresponding con- 
stituents of the phonon modes can be conveniently parametrized by the elementary correlated 
0~ particle-hole operators Aq and Aq defined below. We use the following representation for 
the collective phonon creation and destruction operators A\ oll and A co u 

A \oii = xY / (^) mi c] n c i - yY,(™)i m 4c m , (8) 

mi mi 

and A co u given by the Hermitean conjugate of A\ oll . Here, {0) a b means the matrix ele- 
ments of a single-particle operator O between the nucleon states ip a and ipb- The sym- 
bols a) and a denote the harmonic oscillator "raising" and "lowering" operators a) = 



1 

V2 



{muj) l l 2 r — i{muj) 1 l 2 p 
of the operators 



i 

V2 



[muj) l l 2 r + i{muj) 1 l 2 p 



It is seen that the role 



A o = J2(^)miC m Ci, A Q = ^2(aa) im ,clc m , (9) 

mi mi 

is to create and destroy the correlated particle-hole excitations with the total quantum 
numbers J n = 0~. Indeed, the operator raises the oscillator principal quantum number 
iV by unity, while the role of the spin operator in the internal product is to rotate properly 
the nucleon spins. The combined effect is to remove particles from the orbital with the 
quantum numbers [N, l,j] and place them in the orbital [N + 1,/ ± 1, j]. Thus, due to 
this construction the operator A^ can only create the particle-hole pairs of 0~ type (m and 
i denote the states above and below the Fermi level, correspondingly). The effect of the 
annihilation operator Aq is just the opposite to that of Aq. 

The expressions for the phonon creation and destruction operators A* and A are an 
ansatz for the more general expression of collective RPA phonon operators given in Eq.(^). 
The form of the ansatz is based on the simple commutator properties for the operators (ar) 
and (ap) [|J. with the interactions Vq and V\. (Note that if we substitute the interactions 
Vq and V\ by factorized expressions Vq — > const(ar) a b(ar) d, and Vq — > const(ap) a b(&p) c d, 
the ansatz (H) will be an exact solution of the RPA equations (||)). 



The use of Eq. (||) allows us to get a closed system of RPA equations for the decoupled 
collective phonon modes (T = 0,1), and instead of system of a large number of coupled 
equations (||) we obtain the simplified, approximate sytem of the equations for the amplitudes 
x and y 



EcouAlu ~ [Al n , H ] + ([Al u , V + K]), 



E sp A\ ~ [it H 



(10) 



From the first one of Eqs.flTUD, the closed equations for the collective phonon amplitudes x 
and y can be easily obtained using the following relations 

Po w, 



([Al,V ] 



\U n \ 2 



(A + Al) 



([AlV 1 ]) = -2mpuj(Al-A c 



(11) 



where Vo and V% denote the interaction terms {&i&2)8{r\ — ^2) an d (<Ti(?2){pi, {p*2, <5(r*i — ^2)}} 
respectively expressed in the second quantization form (the curly brackets denote an anti- 
commutator). Hermitean conjugation gives the corresponding relations for the commutators 
of Aq operators. In the first one of Eqs.flTTD we assumed the proportionality of the nucleon 



density p and the potential U: p ~ — u^U, where p and U are the values of the density 
and the potential in the center of the nucleus. 

All the remaining phonon degrees of freedom, marked by the subscript 'sp' in Eqs. (|1C|) , 
are assumed here to be essentially noncollective. The contribution from the interaction 
term in the second one of Eqs.(|T0"D is dropped. The frequencies of those predominantly 
single-particle modes are not shifted from the shell-model value u. 

For the case of two kinds of nucleons (p and n), the collective phonons operators are 
sought in the form A\ oll = x p Aq(p) + x n Ao(n) —y p Ao(p) — y n Ao(n) instead of (H). From (|10|), 
we obtain the following system of equations for the amplitudes x p ,x n , y p and y n 



( 



Ex 






Ex 



CpmZ 



V 



, . CpmZ (-) 



CpmN J,-) 
CpmN n (+) 



— Ex 




u \u\aH 



\ 



-(l+^ (+) ) 



— Ex 



( 



■Ep Up 

•Ep Vp 

•En Vn 



\ 



0. 



where g {±) = g ± g' and g[ ±] = g x ± g[ 



It is convenient to analyze in detail the case of equal numbers of protons and neutrons 
(N = Z). From the consistency condition for (|T2| ) we find the two eigenfrequencies of the 
collective modes Et=o and Et=i to be 



Et=o — w 



Cpm 



Cp XI 1 / 



2 -9i 1 + Tm 9 
Pf > V \U\ 



Et=i — oj 



Cpm A ( Cp , 

1 + I 1 + 



1/2 



, (13) 



The eigenvectors of the system (|T^) are easily found using Eq.fJfJ]), in a form of phonon 
creation operators 



At 



T=0 



CO 



E T=0 (1 + ^ 9l )K 



p'f 



coll 



1/2 



E- 



T=0 



+ 1 + 



CO 

Et=o 



Cpm 
Pf 



(At(p)+Al(n)) 



co 



Cpm 
+ 1 + -Wi 



P 2 f 



(i (p) +A (n)) 



- - 



CO 



nl/2 



co£Z 



Et=1 



+ 1 + 



Cpm , 



2 "fl'l 



- 4(n)) 



+ 



Et=i 



CO 



+ 1 + 



Cpm / 



(Ao(p)-Ao(n)) 



(14) 



where A/" co « is the colective phonon normalization constant. The operators of the independent 
RPA modes should obey the orthogonality and normalization conditions, thus 



[4, A%,] = 0, (0 + [A T , ATr,]0 + ) = 5 TiT ,. 



t in+\ 



(15) 



where the expectation value is taken over the uncorrelated ground state. The first set of 
Eqs.fllBl) (orthogonality) holds by construction, while the second one determines the nor- 
malization. Substituting ([14] ) into ( |i"5|) gives us the normalization relation 

(o + |[i T ,4,]|o + ) = 5 T , T ,(R P + RnWJi- 

Here, for a given kind of nucleons (p or n), the quantity R is 

R = ([(5a), (*#)]) = (£(3 + 2(a?) ii )c\c i ) = 3J\f 



where I is the single-particle orbital angular momentum. The expectation value here is 
reduced to the summation over the occupied states % and gives three times the number of 
particles of given kind, N '. Note that for the spin saturated shell-model ground state of a 



doubly magic nucleus, the term proportional to (al)u does not contribute. Therefore, the 
phonon normalization constant (for N = Z = A/2) is 

Koii = 3 A (16) 

The matrix elements of the operators (af) Ptn and (ap) PjTl between the correlated ground 
state |0 + ) and the excited single-phonon states |0~,T = 0,1) = Ay. =01 |0 + ) can be easily 
obtained by expressing the operators in terms of the phonon creation and destruction op- 
erators Ap =0 1 and At=o,i [see Eqs. ©,(0)]. The results for the matrix elements of the 
isoscalar W P = (0 + |((xp)i|0~, T = 0),W PT = (0 + |(or)i|0~, T = 0) and isovector parity 
violating operators W[ = (0 + \(ap)l\0~ ,T = 1), Wf T = (0 + 1 (c?r) 1 10~, T = 1) are found to 



be: 



wr = r— T 1 , w n p Y mET=oR ^ 



"V 2mE T=0 J ' rr ° ~ \2(l + ^ 9l 



Pp 

To reveal the collective enhancement factors, we compare these results to the matrix 
elements of the corresponding operators between shell-model states. For the (Bp) and a 
given lparticle-lhole state, cJ n Ci|0 + ), we have the following relations 

<0+l(5fl 4,«,|0+> = -i (^f (3SU ~ (^f yfT, ~ {™f W (18) 

where the last expressions indicate the order of magnitude. Np denotes the principal oscil- 
lator quantum number at the Fermi energy of a doubly magic nucleus. We used in (IE) the 
fact that Np ~ (3A) 1 / 3 . Similarly: 

1 \l/2 / i n1/2 / 1 a/2 



4W = ( S ) (-),- ~(^) V^-(^) (^)" 6 da) 

Comparing the matrix elements ( |T7|) to the single-particle ones ([18]), (0) and taking into 
account Eq.( ([T3|) , we see that the matrix elements of the P-odd potentials connecting the 
correlated ground state to the single-phonon states are enhanced by the factors 
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P |(o+|(gp)|o-,r^o)| / w c P 1/3 



and 



^ J(oW.r = o )L + ^ (SA)V> _ (21) 

V A) 



Analogous expressions can be written for the T = 1 matrix elements by substituting g'(g[) 
instead of g{g x ) in (|0|j2l|). 



The total collective enhancement factor F co u has a natural form typical for the low-energy 
collective modes. On the right hand side of Eqs.(pO|),(pTD, the second factor is the collective 



"phase volume" enhancement, F pv , which is roughly equal to the square root of the number 
of single-particle constituents of the phonon, i.e., approximately the number of nucleons in 
the upper major shell, Af), 

F pv ~ Ja~ f ~ N F ~ (3A) 1/3 . 



The first factor in Eqs.(PDD,(PT|) is the "adiabaticity parameter" It reflects the 

change of the frequency of the collective motion compared to the frequency of single-particle 
motion u 



ad 



Note that for the low-lying CT-mode (T = 0), F aa < > 1 gives additional enhancement. 

We turn now to a quantitative discussion of the above results. Using the conventional 
constants of the Landau-Migdal interaction (0),© g = 0.575, g' = 0.725, g\ = —0.5 and 
g[ = -0.26 (see e.g. Refs. 0], p|, |TJ, [pj), we find from Eq.(0): 



E T=0 = 1.03a;, E T=1 = 1.20a;. (22) 

The dependence of Ex on the total strength of the interaction constant is shown on Fig.l. 
It is seen that due to the negative value of g±, g[, the effect of the velocity dependent 
interaction is to reduce the frequency of the isoscalar mode pushing it towards instability 



||. However, the velocity independent component acts in the opposite direction. The 
excitation energy of the isovector mode that is governed by the primed constants tends 
to increase with the interaction strength. (These tendencies are in agreement with those 
found in numerical studies || and with results for the meson exchange interactions Q). The 
combined and opposite action of Vq and V\ on tends to moderate energy shift from u, 
so the values (22) do not mean that the correlations in the ground state are small. The 



collective enhancement of the matrix elements fll7D are F p ~ 10 — 15 for the above values 
of g. 

We should stress that it is not the prime purpose of this work to achieve agreement with 
experiment. The mean field that provides the basis of the single-particle energies is not 
treated in this work. Therefore effects of consistency between the p — h interaction and the 



mean field are not taken into account. Also, we do not include the one-body spin-orbit 
interaction which for 0~ mode could be important especially in lighter nuclei 0. 

As we mentioned, the velocity dependent Landau-Migdal interaction can be viewed as 
a model of a more realistic momentum- dependent n- and p-meson exchange interactions. 
Using V n+P of Eq.(|5|) as a residual interaction within the above formalism we find 



-Et=o — ^ 
Et=i — oj 



fc\A Cp , 



1/2 



1/2 



(23) 



Here, the constants k and <fi and £ are related to the meson-nucleon coupling strenghts: 
k = 2p mq results from the evaluation of the exchange terms with q = 67r(^ r H /7 7r — 

A f 2 

f^rWp) and the nonlocality factors W (W^p — > 1 for m njP — > oo) are W n = 0.11, W p = 
0.69 for the pion and p-meson correspondingly ||. results from the direct terms, = 
— v = -tt^tTO 2 (4 + while the term £ is due to the exchange terms and 

£ = ^.M-l-A-W'^ + 2&W' P ) where W' n = 0.56 and W' p = 0.71. The pion constant is 



equal /J = 0.08, while for the p-meson the coupling constant is in range from j 2 p = 1.86 
("weak p-meson") to f 2 = 4.86 ("strong p-meson") JTSJ, p7[ . 

The phonon operators Aj,, Aj- are given now by the same expressions as before [Eq.flUDl 
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c^ gx ! _ 3k/2 and J + ^ 
given by (]23|) . Correspondingly, the matrix elements of the symmetry violating operators 



but with the replacements 1 + ^™g x -> 1 - 3A;/2 and 1 + &5rd -> 1 + fc/2 and with £ r=0 i 



PT _/ (l-3fc/2)iA V 2 P _ / m£ T ^ y/2 
^° "I 2m£ T=0 J ' W ° ~ \2(l-3k/2)) ' 

^ -l2mE T=1 J ' ^ " \2(1 + A;/2)J ' (24) 

The results for the frequencies are plotted in Fig.l. It is seen that the effect of the 
meson exchange interaction is to push down the isoscalar frequency Et=o, mainly due to the 
7r-meson contribution. At the same time it shift upwards the isovector frequency. Similar 
behaviour of the Et=i was obtained in a more realistic RPA calculations in Ref. ||. Fig.2a 
shows Et as a function of j 2 p . Instability of the 0~,T = occurs for the "weak p-meson", 
as was discussed in [|J using general properties of the linear response |Tj|]. It is interesting 
that the collapse occurs also for the "strong" p-meson, while the stable energy range E ~ 
(0.2 — 0A)u> occurs for 1.86 < f 2 < 4.86. Fig. 2b presents the enhancement factors F p and 
F PT of the weak matrix elements modified according to (p4[) . The collective enhancement 
for the P-odd matrix elements rises from 10 — 20 to considerably larger values as E^ =0 — > 0. 

The enhancement of the P-odd matrix elements connecting the ground state to the single- 
phonon states found here can result in observable P-effects [e.g., in the electromagnetic 
transitions (multipole mixing) |2(J]. For the P-odd (T = 0) matrix element, where the 
enhancement is maximum, the quantity W p oll can reach, for example for 206 Pb, a value of 

W£i = (0 + \g™^^{(a p p p ),p} + g W-^{(a n p n ),p}\0-,T = 0) ~ 20 - 30eV, 

for the realistic weak interaction strengths g p ~ 4 and g n ~ 1 (see, e.g., ||19|| ) consistent 
with the standard parameters of the two-body weak interaction ||21|| , f2"0|| (G is the Fermi 
constant). 

The collective of 0~ states has been already considered in the doorway mechanism ap- 
proach H in the problem of P-odd mixing in compound nuclei. The enhancement in the 
parity violating matrix element due to the collectivity of the 0~ was one of the assumptions 
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in the above works. The results of the present work confirm their idea and the collectivity 
of the P-odd matrix element is a natural result of the RPA treatment of the 0~ states. 

One of us (N. A.) is grateful to B. Mottelson and L. Zamick for discussions. O.V. is 
grateful to V. G. Zelevinsky for discussions of various problems of collective states. The au- 
thors acknowledge the participation at the ECT Trento workshop. The work was supported 
by the US-Israel Binational Science Foundation and the Israeli Science Foundation. 
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Figure captions 



Fig. 1. The 0~ energies Et=o and Et=i [Eq.(|T^), in units of u] as a function of the 
strength of the Landau-Migdal interaction C [taken in units of the standard value (Eq.|3|,|])] 
(solid curves). The dashed curves show Et=q and Et=i as functions of the pion coupling 
constant f% (the scale is given in the upper line of the box) for the p-meson constant j 2 p 
taken in the "weak meson" limit. The collapse of the RPA frequency Et=q occurs when the 
actual value of the pion constant f% « 0.08 is used. 

Fig. 2. a) Energies of the 0~ states E T=0 and E T=1 [Eq.(E3|)], in units of u are plotted 
as a function of the p-meson strength j 2 p calculated for the standard value of the 7r-meson 
constant. Collapse points in E T=0 occur at values of f p approximately equal to the "strong" 
and "weak" p-meson couplings. 

b) The enhancement factors F p (crosses) and F PT (solid lines) [Eq.(p0|)1 as function of 
fp (for the standard value of /|). 
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